Complexity of thermal states in quantum spin chains 
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We study quantum correlations and complexity of simulation, characterized by quantum mutual 
information and entanglement entropy in operator space respectively, for thermal states in critical, 
non-critical and quantum chaotic spin chains. A simple general relation between the two quantities 
is proposed. We show that in all cases mutual information and entanglement entropy saturate with 
the system size, whereas as a function of the inverse temperature, we find logarithmic divergences 
for critical cases and uniform bounds in non-critical cases. A simple efficient quasi-exact method 
for computation of arbitrary entropy related quantities in thermalized XY spin chains is proposed. 

PACS numbers: 02.30.Ik, 05.70.Fh, 75.10.Pq, 03.67.Mii 



I. INTRODUCTION 

Physical properties of non-integrable strongly corre- 
lated quantum systems are notoriously difHcult to calcu- 
late due to exponential growth of Hilbert space dimension 
with the system size. Nevertheless, many numerical and 
analytical approaches have been developed that provide 
results for special properties and specific systems. The 
question arises, which quantities are amenable to classi- 
cal calculation in spite of generic system complexity? 

Recently, two results emerging from the field of quan- 
tum information shed light on the question. First, it 
has been shown that for one-dimensional (ID) systems 
the entanglement content of ground states grows loga- 
rithmically with the number of particles for fermionic 
critical systems and saturates for non-critical or bosonic 
systems |1„ 2,, j5„ . Second, density matrix renormaliza- 
tion group (DMRG) method ^ has been reinterpreted 
and optimized within matrix product state (MPS) ansatz 
@) 0] for the many body wave-function. Its complex- 
ity is essentially given by the entanglement entropy of 
the system's bipartition. Despite being essentially the 
best method available it is not known precisely in which 
instances is DMRG computationally efficient, for exam- 
ple it has been shown to be inefficient for computing 
the evolution of homogeneous ID non-integrable systems 
in real time Q but becomes efficient in the presence of 
disorder due to the emergence of (many body) localiza- 
tion Thermal states at finite temperature have been 
simulated using an extension of MPS to density opera- 
tors {matrix product operators MPO), and by performing 
the evolution in imaginary time (inverse temperature) 
Q- However, no analysis of efficiency scaling with the 
system size n and inverse temperature /3 have been per- 
formed. Note that finding a ground state - corresponding 
to the limit /3 — > oo - is proven to be a hard problem [l^l , 
even though it can often be performed efficiently Q . For 
a different algorithm for computing thermal states see 
Ref. pjj. In Ref. [l3| an upper bound on the necessary 
MPO dimension is found for thermal states which un- 
fortunately scales exponentially with (3. It is not clear if 
one can do better, that is, if finite temperature states can 
be calculated efficiently, uniformly in /3? In the present 



work we show that this is indeed possible. 

Complexity will be characterized by the operator space 
entanglement entropy (OSEE) of thermal states 

which is directly linked to the minimal necessary rank 
of MPO ansatz. OSEE is defined as von Neumann en- 
tropy of the reduced density matrix defined with respect 
to a physical density operator treated as an element of a 
Hilbert (Fock) space of operators. In addition we shall 
show that OSEE behaves in the same way as the quantum 
mutual information (QMI) Ia-.b 
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which is a natural 
measure of quantum correlations in a bipartition A + B 
and an upper bound for distillable entanglement. Our 
main result is that OSEE and Ia-.b in thermodynamic 
limit n — > OO grow as ~ (c/3) log2 (3 for critical (gap- 
less) systems and saturate for non-critical (gapped) sys- 
tems, while they always saturate with f3 for a finite n. 
This shows that the linear in f3 upper bound on QMI 
of Ref. [T^ is far from optimal and DMRG computations 
of large (infinite) critical ID quantum systems should 
be efficient at all nonzero temperatures. This has to be 
contrasted with typical linear growth of OSEE in t in 
real time evolution Q. Note that OSEE and QMI al- 
ways saturate with n, for any fixed and hence scale 
very differently from thermodynamic entropy of a block 
of length n, which for a critical system behaves as [l^ 
S" (X ^ for large n. Our results are based on analytical 
and numerical calculations in several examples of inte- 
grablc and non-integrable spin chains. We outline a novel 
method for quasi-exact calculations of OSEE and Ia-.b in 
XY spin 1/2 chains using the concept of Fock space of op- 
erators [131, whereas we rely on numerical MPO-DMRG 
simulations for the other cases. 



II. METHODS 

First, we describe how to construct a thermal state 
and compute OSEE in a general spin chain, re- 
gardless of its integrability, using MPO ansatz and 
imaginary time evolution. A thermal state p((3) = 
exp {—j3H) I tr [exp (— )] of Hamiltonian if of a chain 
of n spins 1/2 can be expressed as an element of 4" 
dimensional Hilbert space /C, \p) = ^gCs\cr-) spanned 



2 



by Pauli operators a- = (tJ^---ct^", s = si,...,s„, 

S„, e {0,1,2,3}, Cr° = l,Crl'2^3 = j^^^^^ pj,^^, 

uct in K is defined as {x\y) = 2^"trx'''2/. Within MPO 
ansatz expansion coefficients Cs are represented in terms 



of in D X D matrices Af 



1, 



tr (Af 



0,1,2,3, as 
(1) 



For exact MPO representation matrix dimension of A^ 
must be equal to the number of nonzero Schmidt coef- 
ficients for a bipartite cut of a vector Cg at the corre- 
sponding site. Although for generic density operators p 
the required dimension D is exponentially large (in n) in 
practice one truncates the matrices to a smaller dimen- 
sion D and uses ([1]) as a variational ansatz. As a rough 
estimate of a minimal D ^ 2^ , which is a direct measure 
of complexity of p, one can use OSEE S*' defined as an 
entanglement entropy of a super- ket \p) 

5« = -tr^(Rlog2R), R={p\p)-Hi-b\p){p\ (2) 

where a subscript A denotes the first ha sites and B its 
complement. More detailed analysis [l^ specifies exact 
criteria for the applicability of matrix product ansatz ([T|) 
in terms of properties of the moments of Schmidt coeffi- 
cients. Writing the expansion coefficients as a 2"-* x 
2n-riA dimensional matrix Cs s = Cc s , the reduced 
density matrix ^ is given by R = CC^/tr (CC^) and 
quantifies the bipartite correlations and the difficulty of 
a MPO simulation of a given density operator p. For 
use of OSEE in a non-MPO context see Ref. [H. S"" is 
directly related to the entanglement entropy of physical 
states only in the case of pure states, p = ji/') , in which 
case 5" is just twice the entanglement entropy of a pure 
state ip, e.g. the ground state at zero temperature. How- 
ever, even though 5" is not an entanglement monotone 
for general p it can be used as a bipartite entanglement 
estimator for thermal states p{/3). Namely, as we will 
demonstrate later, its scaling on n and (3 is the same as 
that of QMI Ia:B, defined as 



Ia:B = S{pa) + S{pb) - S{p), 



(3) 



where pA = tr^p, Pb = tr^ip, and S{a-) ~ — tr (cr log2 cr). 
QMI measures total correlations between A and B and 
is an upper bound for relative entropy of entangle- 
ment which in turn is an upper bound for distillable 
entanglement Ed{p) [2^. MPO representation for 
a thermal state apart from normalization, is 

computed by performing imaginary time evolution in 
/?, p{l3 + e) = exp(-ei?/2)p(/?)exp(-ei?/2), starting 
at an infinite temperature state p{0) = 1/2" having 
D = 1. Such /3-evolution increases OSEE and as a 
consequence, dimension D has to increase as well in 
order to preserve the accuracy. Assuming the empirical 
fact that a fixed small e (in our case e = 0.05) is enough 
to produce accurate results in a finite number of steps 
(for all (3 < oo), we conclude that time complexity of 
simulation is polynomial in D and /3. Note that in 



contrast to a (unitary) real time evolution, imaginary 
time evolution does not preserve canonical Schmidt 
decomposition structure of matrices A^" In or- 

der to improve stability we apply local rotations every 
few steps in order to "reorthogonalize" the matrices A^™ . 

Quasi-exact calculation of and Ia-.b for thermal 
states is possible for intcgrable chains which can be 
mapped to quadratic fermionic models by Wigner- 
Jordan transformation, e.g. ID XY model. Most 
general such situation is described by a Hamilto- 
nian H ~ w ■ Hw = '}2j"i=i where H is 

an anti- symmetric Hermitian 2n x 2n matrix and 



W2m-1 



n 



W2m = CTmUZi^ '^i are Majo- 
rana operators satisfying Iwj , w/| = wjwi +wiWj = 25 ji. 
It is crucial to recognize [IJ, [l3| that K. is in fact a Fock 
space if spanned by 4" products of wj , Pa = w^^ ■ ■ ■ ^2??" j 
aj G {0, 1}, and equipped with adjoint Fermi maps 
Cj\Pa) = Saj^i\wjPa). In our notation • designates a 
linear map over the Fock space of operators /C. Note 
that OSEE is just an entanglement entropy of an 
element of K.. Therefore, using the formalism of 0] lifted 
to an operator space, one can calculate S'" from the 
anti- symmetric correlation matrix T defined by 



1=1 



(e-'3g|apa,|e-^^) 

lp-PH\p-f3H\ 



= S 



PI 



■iP 



p,g=l,...,4n (4) 



where 02^-1 = {cj + c]), a2j = i(cj - c]), j 
satisfying {ap,aq} 



1. 



,2n 



written in terms of dynamics 



'pq. Thermal state in ([4]) can be 

(5) 



using a quadratic adjoint Hamiltonian map Ti. over /C, 
n = (c + c') ■ll{c + c'), since (cj +c])|p) = \wjp). Di- 
agonalizing the matrix H we obtain eigenvalues A^, —Xk, 
ordered such that < Ai < . . . < A„, and the correspond- 
ing orthonormal eigenvectors w^,, in terms of which we 
write the adjoint Hamiltonian in ([5]) in the normal form 
^ = ELi 4Afc(6l.6fc - i) where bk = vl ■ (c + c^)/V2. 
Identities (l|6fe|l) = 0, {l\bkbi\l) = 0, (l|6^6i|l) = ^Su 
allow us to employ the Wick theorem to evaluate ([4]) 



r2j-i,2/-i = -r2j,2/ = -2y^ tanh(4/?Afc)Im(ffc ;) 

k=l 



r2j,2/-i = Sji - 2^(1 



1 



cosh(4/3Afc) 



)Mvl,v^,i)- (6) 



Finally, OSEE is obtained from the eigenvalues iii^j of 
An A X Aha upper-left block of T as = Xlj"! H2{^) 
where H2{x) = — xlogj a; — (1 — a;) log2(l — x) 0. 

Similarly, we can also easily compute the block en- 
tropy S{pa) of a thermal state p(f3) by diagonalizing 
H = X]fc=i 4Afe(6|6fc - i) with operators bk = vl ■ w/ V2 
(which are elements of /C and not maps over /C like 
bk) and computing the correlation matrix F', defined as 
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tr [wjWie-l^")/ tve-l^" = Sji + ir^.,, j,l = l,..., 2n. Ex- 
pressing Wj and H in terms of b^, bj,, and using the Wick 

theorem with tr(6[.&i) = ^6ki trl results in a general 
2-point thermal correlation function 



r' 



k=l 



(7) 



Diagonalizing 2nA x 2nA block of T' corresponding to 
the sub-lattice A, yielding eigenvalues ±ii^j, we obtain 

the block entropy S{pa) = J2]=i -f^2(— j-^)- Considering 
also B, and A + B instead of A we get QMI 



III. RESULTS 

We shall show results for OSEE and QMI for several 
families of critical and non-critical as well as integrable 
and non-integrable quantum spin chains 



HxY = 2^ [ -^^i ^i+i + 



1=1 

n-1 



1=1 1=1 

ri-l 



1=1 

(8) 

n 



I I 



(=1 



1=1 



all for open boundary conditions, and a symmetric cut 
at UA = n/2 where largest bipartite entanglement is ex- 
pected. Majorana representation of the integrable XY 
case which is necessary for explicit calculations [Eqs. jH- 

I?])] is HxY = -^J2?=l {^W2lW2l + l- ^W2l-lW2l+2) - 
^J2?=l hW2l-lW2l- 

In Fig. [T] we show the results for quantum critical sys- 
tems for which the entanglement entropy of the ground 
state grows as So ^ (c/3) logj n Q where c is the central 
charge of the corresponding conformal field theory 
However, we take several very large n and plot f3 depen- 
dence of OSEE and QMI, for (a) XX model Hxril = 
0, /i = 0), (b) transverse Ising model Hxy{i = 1, /i = 1), 
and (c) XXZ model Hxxzi^ = 0-5, hi ee 0). These 
critical models are all completely integrable, however 
we were able to make exact calculations only for (a,b), 
whereas we used MPO-DMRG simulations for (c) where - 
since evaluating QMI ([3|) is computationally difficult even 
within MPO ansatz - we have instead computed mutual 
purity, defined as Ip = \og2 [P{p)/{P{pa)P{pb))] with 
P((t) = tr cr^. It is well known that purity often well 
describes the behavior of von Neumann entropy and its 
complement is sometimes termed linear entropy. In addi- 
tion we checked for thermal states of small systems that 
mutual purity and QMI behave similarly. 

From Fig.[l]we conclude that in the large n limit OSEE 
grows logarithmically with /? 



S\p)^-l0g2p + c', 



with parameters c = 1 and c' = 1.17 for XX model (a), 
c = 1/2 and c' = 0.75 for transverse Ising model (b), and 
c — 1 and c' — 1.63 for XXZ model (c). c is exactly the 
central charge of the corresponding model. Presently we 
do not have an explanation for this finding and we state 
it as a conjecture. For fixed n, S*" saturates for (3 larger 
than the inverse spectral gap /3* = 1/AE at twice the 
ground state value 25*0 ~ (2c/3) logj n. 
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FIG. 1: Logarithmic growth of OSEE S" (upper curves) and 
QMI Ia:B (lower curves) for thermal states of critical models: 
(a) XX chain Rxviy = 0,h = 0), (b) transverse Ising chain 
Hxri^y = = 1), both for n = 1024,256, and (c) XXZ 
chain Hxxz{A = 0.5, hi = 0) for n = 256,32,8. In case 
(c) we show mutual purity Ip instead of Ia.b- Dotted lines 
indicate conjectures ([9]), and pO|) combined with ^ [see text]. 

Next, let us have a look at non-critical systems. First, 
we will chose two non-integrable spin models that dis- 
play typical signatures of quantum chaos (Fig. (a) 
Heisenberg XXZ model in a staggered magnetic field 
Hxxz{A = 0.5, hi = -(1 + (-1)0/2), and (b) Ising 
model in a tilted magnetic field Hi{hx — l,hz = 1). We 
can see that for non-critical systems which possess a fi- 
nite energy gap AE between the ground and first excited 
states, OSEE saturates at n and (3 independent values for 
sufficiently large /3 ^ /?*, so the classical simulation of 
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FIG. 2: Saturation of OSEE S" (upper curves) and mu- 
tual purity Ip (lower curves) for quantum chaotic noncnt- 
ical systems, namely (a) for XXZ model in a staggered field 
Hxxz{A = 0.5, hi = -(1 + (-l)')/2) and (b) for tilted Ising 
Hi{h^ = = 1), for n = 8,32, 128. Note that the curves 
for n = 32 and n = 128 are practically indistinguishable. 
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FIG. 3: OSEE 5" (full curves) and QMI Ia-.b (broken curves) 
for noncntical XY chain with 7 = 0.5, h = 0.9, and n = 
20, 40, 60. Note a crossover between two saturation plateaus 
due to a single exponentially close (in n) excited state. The 
inset shows h dependence of the gap decay length ^, computed 
from the fit to A E (x e~ "/^, torn = 30, . . . , 60. Vertical line 
indicates h* = a/1 — 7^. 

quantum thermal states for gapped systems is even more 
efficient than for critical ones. Although this might seem 
expected, one should know that for OSEE we do not have 
theoretical uniform bounds in n like for QMI [lH|- Note 
that in the presence of weakly broken symmetries, the 
gap AE between two lowest states can get exponentially 
small with n even though the system is not critical. This 



for instance happens in XY model. In Fig. [3] we show 
results for a non-critical integrable XY model ([8]) with 
7 = 0.5 and h = 0.9, where /3* = 1/AE grows exponen- 
tially with n. Actually we have shown for any Hxy, that 
AE = 4Ai where Ai is the smallest eigenvalue of matrix 
H ([5]), and which in turn has been found to decay expo- 
nentially Ai ^ exp(— n/^). Interestingly, for /i^ +7^ < 1, 
the gap decay length ^ is practically insensitive to pa- 
rameters 7, h, whereas at the point h* ~ \/l — 7^, where 
the ground state is doubly degenerate and separable [2l[ , 
^ instantly starts increasing with increasing h—h* (inset 
of Fig. [3]). However, there is no phase transition at h* . 

In all cases studied QMI behaves in essentially the 
same way as OSEE. In fact, for models solvable by 
Wigner- Jordan transformation Ia-.b {(3) almost overlaps 
with evaluated at /?/4 

Ia:bW) ~ S^P/A) (10) 

apart from finite size effects at small and large /3 
(Figs. Ill3p . We conjecture that pI7|) becomes asymp- 
totically exact for large n. Factor 4 may be understood 
from comparing the formalisms of Majorana operators 
Wj and the adjoint maps Cj. Since QMI is an upper 
bound for distillable entanglement [2^ this means that 
exact quantum bipartite entanglement for thermal states 
is also small, i.e., it is either upper bounded by n and 
/?-independent value for non-critical systems (Figs. I2I3P 
or it is upper bounded by n-independent value ~ log2 (3 
for critical systems. 



IV. CONCLUSIONS 

By studying several typical models, we have demon- 
strated that classical simulations of thermal states of 
interacting quantum systems in ID are efficient. Us- 
ing DMRG simulations combined with exact calculations 
for XY models we found that in non-critical systems 
complexity and bipartite entanglement measures are uni- 
formly bounded, both in system size n and inverse tem- 
perature /3, whereas they exhibit a universal logarithmic 
divergence in inverse temperature for quantum critical 
systems. 

Our results represent a uniform finite temperature ex- 
tension of the ground state area-law 

^,U,m- Specif- 
ically, for critical systems, logarithmic divergence of 
ground state block entropy with the size of the block 
is reflected in logarithmic divergence in inverse tempera- 
ture (3 of the bipartite entanglement of thermal (mixed) 
states in the thermodynamic limit. Even though our re- 
sults are not rigorous they suggest optimal bounds on 
entanglement and complexity of thermal states that are 
qualitatively stronger than previously known (Tsj . 
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